REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  the  burden,  to  the  Department  of  Defense,  Executive  Services  and  Communications  Directorate  (0704-0188).  Respondents  should  be  aware 
that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  0MB 
control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION. 


1.  HEPOKX  DAJE  (DD-MM-YYYY)  2.  REPORT  TYPE 

14-08-2014  Journal  Article 


4.  TITLE  AND  SUBTITLE 

Inserting  Tides  and  Topographic  Wave  Drag  into  High-resolution  Eddying 
Simulations 


3.  DATES  COVERED  (From  -  To) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Brian  Arbic,  Maarten  Buijsman,  Eric  Chasignet,  S.T.  Gamer,  John  A  Goff, 
S.R.  Jayne,  E.  Joseph  Metzger,  James  G  Richman,  Jay  Shriver,  Patrick  G 
Timko,  D.S.  Trossman,  Alan  J  Wallcraft 


0602782N 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory 

Oceanography  Division 

Stennis  Space  Center,  MS  39529-5004 


73-8677-03-5 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

NRL/PP/7320-14-21703 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 

One  Liberty  Center 

875  North  Randolph  Street,  Suite  1425 

Arlington,  VA  22203-1995 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  is  unlimited. 


14.  ABSTRACT 

Tidal  forcing  has  been  inserted  into  eddy-resolving  global  simulations  of  the  HYbrid 
Coordinate  Ocean  Model  (HYCOM).  The  tidal  flows  over  realistic  topography  generate  internal 
tides^  i.e.  internal  waves  at  tidal  frequencies.  Similarly^  geostrophic  flows  over  rough 
topography  generate  internal  lee  waves.  While  it  is  not  yet  possible  to  resolve  internal  lee 
waves  in  global  models^  their  impact  can  be  parameterized.  We  have  inserted  parameterized 
topographic  lee  wave  drag  in  global  simulations  that  do  not  include  tidal  forcing. 


15.  SUBJECT  TERMS 


HYCOM,  tidal  forcing,  topographic  wave  drag,  bottom  drag 


16.  SECURITY  CLASSIFICATION  OF: 
a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 

Unclassified  Unclassified  Unclassified 


17.  LIMITATION  OF 
ABSTRACT 


18.  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 

OF  Joseph  Metzger 

PAGES  - - - r - 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

^  (228)  688-4762 


Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39.18 


PUBLICATION  OR  PRESENTATION  RELEASE  REQUEST 


iV:  E':. ■ 


ill 


Ref:  (a)  NRL  Instruction  5600.2 
(b)  NRL  Instruction  5510.40D 

End:  (1 )  Two  copies  of  subject  paper 
(or  abstract) 


(  )  Abstract  only,  published  ( 

(  ) Book  ( 

(  )  Conference  Proceedings  ( 

(refereed) 

(  )  Invited  speaker  { 

(  )  Journal  article  (refereed)  ( 

(  )  Oral  Presentation,  published  ( 

(  )  Other,  explain 


(  )  Abstract  only,  not  published 


#  ^  - 


(  )  Book  chapter 

(  )  Conference  Proceedings 

(not  refereed) 

{  )  Multimedia  report 

( X  )  Journal  article  (not  refereed) 

(  )  Oral  Presentation,  not  published 


STRN  NRUJA/7320-14-2170 

Route  Sheet  No.  7320/ _ 

Job  Order  No.  73-  8677-  03-  5 

Classification  X  U 

Sponsor  ONR  _ 

approval  obtained  _  yes 


Title  of  Paper  or  Presentation 

Inserting  Tides  and  Topographic  Wave  Drag  into  High-  resolution  Eddying  Simulations 

Author(s)  Name(s)  (First, Ml,  Last),  Code,  Affiliation  if  not  NRL 

Brian  ArbiC  Univ.  of  Michigan  Maarten  Buijsman  University  of  New  Orleans  Eric  Chassignet  Florida  State  University  ST.  Garner  NOAA  -  New  Jersey 

John  A.  Goff  Univ.  of  Texas  at  Austin  S.R.  Jayne  WHOI  E.  Joseph  Metzger  7323  James  G  Richman  7323  Jay  Shriver7323 
Patrick  G.  Timko  Bangor  Univ.  -  UK  D.S.  Trossman  Univ.  of  Michigan  Alan  J.  Wallcraft  7323 


It  is  intended  to  offer  this  paper  to  the 


(Name  of  Conference) 


/^  /  /  i  'A  ^  Place  and  Classification  of  Conference) 

^ X  cju. /i  ^  ^ f- 

and/or  for  publication  in  Exchanges  --eLIVAR  NewslettcrriJnclassified  _ 

(Name  and  Classification  of  Publication)  (Name  of  Publisher) 

After  presentation  or  publication,  pertinent  publication/presentation  data  will  be  entered  in  the  publications  data  base.  In  accordance 
with  reference  (a). 

It  is  the  opinion  of  the  author  that  the  subject  paper  (is _ )  (is  not _ ^)  classified,  in  accordance  with  reference  (b). 

This  paper  does  not  violate  any  disclosure  of  trade  secrets  or  suggestions  of  outside  individuals  or  concerns  which  have  been 

communicated  to  the  Laboratory  in  confidence.  This  paper  (does _ )  (does  not _ X)  contain  any  militarily  critical  technology. 

This  subject  paper  (has _ )  (has  never  X  )  been  incorporated  in  an  official  NRL  Report. 


E.  Joseph  Metzger,  7323 _ 

Name  and  Code  (Principal  Author) 


Signature) 


CODE 


SIGNATURE 


DATE 


Author(s) 


■Cf2y^£jr 


Section  Head  VI 

. . . 

Branch  Head  ^ 

Gregg  A.  Jacobs,  7320 

Division  Head 

Ruth  H.  Preller,  7300 

Security,  Code 

1231 

Office  of  Counsel, Code 

_ 1008.3 _ 

ADOR/Director  NCST 

E.  R.  Franchl,  7000 

Public  Affairs  (Unclassified/ 
Unlimited  Only),  Code  4 

Division,  Code 

Author,  Code 


HQ-NRL  5511/6  (Rev.  12-98)  (e) 


Need  by 


COMMENTS 


/i  ,T/7/n  /  ^ 


Publicly  accessible  sources  used  for  this  publication 


A  /; 
s/^1m 


6-2')t 


THIS  FORM  CANCELS  AND  SUPERSEDES  ALL  PREVIOUS  VERSIONS 


1 .  Release  of  this  paper  is  approved. 

2.  To  the  best  knowledge  of  this  Division,  the 

subject  matter  of  this  paper  (has _ ) 

(has  never  x  )  been  classified. 

1 .  Paper  or  abstract  was  released. 

2.  A  copy  is  filed  in  this  office. 


pubucahon  or  presentation  rei.ease  request 


hJiSs:, 


R^t  (a)  NRL  fns<ructk>n  5600.2 
(b)  NRt  In^iruotion  S510.40D 

End;  (i )  Two  copies  of  subiect  paper 
(or  abstract) 


)  Abstract  orily,  published 
)Book 

)  Con(erenc<»  Proceedings 
(refereed) 

)  Invited  speaker 
)  Journal  artde  (refereed) 

)  Oral  Presentation,  publi^ed 
)  Other,  explain 


(  )  AhsincX  only,  not  published 

/  )  Book  chapter 

(  )  Conference  Proceedings 

(not  refereed) 

(  )  Multimedia  report 

( X  }  Journal  article  (rtot  refereed) 

(  )  Oral  Presentation,  not  published 


Route  Sheet  fte.  7320/ 
Job  Order  No. 
Cla&sificatkm  A 
Sponsor  ONR  sjj 
approval  obtained 


U 


Title  of  Paper  or  Presentation 

InserUos  Tides  and  Topogi^phlc  WOve  Drag  Into  High'-  resolution  Eddying  Simulations 


Author(s}  Name(s)  (First,Mf,LB^r  Code,  Amatiofi  ifaoit  NRL 

Biiert  AitdcUnlv.ofMidiiBan  Maarten  BuUsmanUw/flraiiyofNewOftean^  Eric  Chassignet  Ftodda  Stats  Universitv  ST.  Gamer  NOAA-NawJerEoy 
John  A.  Goff  Unhr.orTexMatAudffi  S.R  Jayne  WKOi  IS,  Joseph  Metzger  7323  James  G  Riohrnan  7323  Jay  Shriver7333 
Patrick  G.  Timko  Bangor  Univ.  -  UK  D.S.  Trossman  unh'.  of  Micfiigan  Alan  J.  Wallcraft  7323 


U  is  intended  to  offer  this  paper  to  the 


(Name  ofCanfamnca) 


n  /  rvJ  ^  {OK}\a,  Ptece  and  Cfaasl/fco/wn  of  Conference; 

2^ //A-  ^ ^  S'  A/e^CtJx 


and/br  for  publication  in  Exchanges  ■  CyVAR-NewsieitorrJUnciagstfieq  _  _ _ _ _ 

(Na/ne  and  Ctese^jfibn  ofPudUcathri)  "  (Ateme  o/PubiiB/ied 

After  presentation  or  publication,  pertinent  publicatlon/presentatfon  data  will  be  entered  In  the  publfoations  data  base,  In  accordance 
w«lh  reference  (a). 

it  is  the  opinion  of  the  author  that  the  subject  paper  (Is _ )  (is  not  classified,  In  accordance  with  reference  (b); 

This  paper  does  not  violate  any  disclosure  of  trede  secrets  or  suggestions  of  outside  individuals  or  concerns  which  have  been 

communicated  to  the  Laboratory  in  confidence.  This  paper  (does _ )  (does  not _ ^ )  contain  any  militarily  critical  technology. 

This  sub^  paper  (has _ _ )  (has  never  X  )  b<ien  incorporated  in  an  official  NRL  Report 


E.  Joseph  Metzger.  7323 


Name  and  Cods  (Pdnc^a/^lu/hod 


Public  Affairs  (Unciaa$m<U 
Un«mtedOn/y;,  Code  ^ 


DfVleiDn,  Code 


Author,  Cods 

HQ4iHL  S5tl«  (Rev.  IZ-W)  (e) 


THIS  FORM  CANCELS  AND  SUPERSEDES  ALL  PREVIOUS  VERSIONS, 


^0/00  39Vd 


3Pi3Sl3lZ2Z 


PI -13  P12Z/9Z/92 


Exchanges 


No.65  (Vol  19  No.2)  July  2014 


Special  Issue:  High  Resolution  Ocean  Climate  Modelling 


10  km 
resolution 


mesh 

transition 


80  km 
resolution 


CLIVAR  Ocean  &  Climate:  Variability,  Predictability  and  Change  is  the 
World  Climate  Research  Programme’s  (WCRP)  project  on  ocean-atmosphere  interactions. 

WCRP  is  sponsored  by 

the  World  Meteorological  Organization, the  International  Council  for  Science  and 
the  Intergovernmental  Oceanographic  Commission  of  UNESCO. 


WGRP-ii' 


World  Climate  Research  Programme 


Inserting  tides  and 
topographic  wave 
drag  into  high- 
resoiution  eddying 
simuiations 

Brian  K.  Arbid*,  Maarten  C.Buijsman^, 
Eric  P.Chassignet^,  Stephen  T. Garner^, 
Steven  R  Jayne^,  E.  Joseph  Metzger®, 
James  G.  Richman®,  Jay  E.  Shriver®, 
Patrick  G.  Timko^S  David  S.  Trossman^A 
and  Alan  J.  Wallcraft® 

1  Department  of  Earth  and  Environmental 
Sciences,  University  of  Michigan,  Ann  Arbor, 
Michigan,  USA 

2  Department  of  Physics,  University  of  New 
Orleans,  New  Orleans,  Louisiana,  USA 

3  Department  of  Earth,  Ocean,  and  Atmospheric 
Science  and  Center  for  Ocean-Atmospheric 
Prediction  Studies,  Florida  State  University, 
Tallahassee,  Florida,  USA 

4  National  Oceanic  and  Atmospheric 
Administration  (NOAA)/Geophysical  Fluid 
Dynamics  Laboratory,  Princeton,  New  Jersey, 
USA 

5  Physical  Oceanography  Department,  Woods 
Hole  Oceanographic  Institution,  Woods  Hole, 
Massachusetts,  USA 

6  Oceanography  Division,  Naval  Research 
Laboratory  (NRL-SSC),  Stennis  Space  Center, 
Mississippi,  USA 

7  Center  for  Applied  Marine  Sciences,  Marine 
Science  Laboratories,  Bangor  University,  Menai 
Bridge,  Anglesey,  UK 

8  Department  of  Atmospheric  and  Oceanic 
Sciences,  McGill  University,  Montreal,  Canada 

Corresponding  author:  arbic@umich.edu 

1.  Introduction 

We  report  here  on  two  of  our  recent  development  efforts 
for  global  high-resolution  ocean  models,  undertaken  in  the 
HYbrid  Coordinate  Ocean  Model  (HYCOM;  Chassignet  et 
al.,  2009).  The  simulations  discussed  here  are  “eddying” 
(Hecht  and  Hasumi,  2008),  meaning  that  they  include  an 
energetic  mesoscale  eddy  field.  In  one  of  our  efforts  (Arbic 


et  al.,  2010,  2012;  Timko  et  al.  2012,  2013;  Shriver  et  al.  2012, 
2014;  Richman  et  al.,  2012),  we  have  inserted  tides  into 
global  eddying  simulations.  Tidal  flows  in  a  high-resolution 
global  model  with  realistic  rough  topography  generate 
internal  tides-internal  waves  of  tidal  frequency.  Similarly, 
geostrophic  flows  over  rough  topography  generate  internal 
lee  waves  (Bell,  1975).  It  is  not  yet  possible  to  resolve 
internal  lee  waves  in  global  models,  but  the  impact  of  internal 
lee  waves  can  be  parameterised.  In  the  other  effort  reported 
here,  we  have  inserted  parameterised  topographic  lee  wave 
drag  into  global  eddying  simulations  that  do  not  include 
tides,  with  a  preliminary  focus  on  the  resulting  model  energy 
budget  (Trossman  et  al.,  2013). 

Global  simulations  of  the  tides  and  the  atmospherically 
forced  eddying  general  circulation  have  long  been  done 
separately.  Global  general  circulation  models  began  to  resolve 
mesoscale  eddies  in  the  1990’s  (e.g.,  McClean  et  al.,  1997). 
Since  then  eddying  models  have  become  increasingly  realistic 
(e.g.,  Hecht  and  Hasumi,  2008).  Global  modelling  of  internal 
waves  is  a  newer  endeavour.  In  the  first  global  simulations 
of  internal  tides  (Arbic  et  al.,  2004;  Simmons  et  al.,  2004), 
only  tidal  forcing  was  employed,  and  the  stratification  was 
taken  to  be  horizontally  uniform.  The  insertion  of  tides  into 
atmospherically  forced  general  circulation  models  allows  for 
internal  tide  propagation  in  a  realistic  horizontally  varying 
stratification,  for  interactions  between  internal  tides  and 
eddies,  and  for  the  co-existence  of  tides  and  near-inertial 
waves,  another  important  class  of  internal  waves  (e.g., 
Simmons  and  Alford,  2012).  In  global  high-resolution 
models  with  simultaneous  atmospheric  and  tidal  forcing,  the 
internal  wave  spectrum  is  beginning  to  be  resolved,  just  as 
mesoscale  eddies  began  to  be  resolved  in  global  models  two 
decades  ago.  Our  simulations  with  simultaneous  tidal  and 
atmospheric  forcing  are  being  used  to  address  a  variety  of 
scientific  and  operational  questions. 

The  global  energy  budget  has  been  a  topic  of  great  recent 
interest,  largely  because  the  mixing  occurring  when  energy 
is  dissipated  may  exert  a  strong  control  on  the  oceanic 
meridional  overturning  circulation  (Munk  and  Wunsch, 

1998).  The  mechanisms  underlying  the  dissipation  of 
the  eddying  oceanic  general  circulation  are  still  under 
investigation.  Candidate  dissipation  mechanisms  include  the 
transfer  of  energy  from  geostrophic  flows  to  submesoscale 
motions  in  the  upper  ocean  (e.g.,  Capet  et  al.,  2008),  bottom 
boundary  layer  drag  (Sen  et  al.,  2008;  Arbic  et  al.,  2009), 
and  breaking  of  internal  lee  waves  generated  by  geostrophic 
flows  over  rough  topography  (e.g.,  Naveira-Garabato  et  al., 
2004).  Nikurashin  and  Ferrari  (2011)  and  Scott  et  al.  (2011) 
estimated  the  globally  integrated  energy  flux  of  geostrophic 
flows  into  internal  lee  waves  over  rough  topography,  and 
found  it  to  be  a  substantial  fraction  of  the  ~1  TW  of  energy 
put  by  the  wind  into  geostrophic  flows  (e.g.,  Wunsch,  1998, 
and  others).  Both  the  Nikurashin  and  Ferrari  (2011)  and 
Scott  et  al.  (2011)  estimates  were  “offline”;  they  utilised 
bottom  flows  from  eddying  simulations  that  did  not  employ 
topographic  wave  drag  as  they  ran.  In  Trossman  et  al.  (2013) 
we  inserted  inline  topographic  lee  wave  drag  into  eddying 
HYCOM  simulations,  thus  ensuring  feedbacks  between 
the  model  bottom  stratification  and  flow  fields  and  the 
topographic  internal  lee  wave  drag. 

The  HYCOM  simulations  shown  here  are  forced  by  output 
from  the  Navy  Operational  Global  Atmospheric  Prediction 
System  (NOGAPS);  going  forward  the  simulations  will  be 
forced  by  the  Navy  Global  Environmental  Model  (NAVGEM). 
The  simulations  shown  here  are  run  on  a  combination 
Mercator  and  tripolar  grid  with  an  equatorial  grid  resolution 
of  0.08°  (8.9  km)  and  are  coupled  to  a  sea  ice  model  at  high 
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latitudes.  Global  HYCOM  simulations  at  0.04°  (4.5  km)  have 
also  been  performed,  but  are  not  discussed  here  (see  article 
by  Chassignet  et  al.  2014,  this  issue).  The  model  reproduces 
the  general  circulation  of  the  global  ocean,  the  strength  and 
variability  of  the  western  boundary  currents,  such  as  the  Gulf 
Stream  and  Kuroshio,  and  the  mesoscale  eddies  generated 
by  instabilities  of  the  major  currents  (Thoppil  et  al.,  2011). 


2.  Insertion  of  tides  into  eddying  simulations 

In  our  HYCOM  simulations  forced  by  both  atmospheric 
fields  and  tides  (Arbic  et  al.,  2010,  2012),  we  include  the 
four  largest  semidiurnal  tidal  constituents  (M2,  S2,  N2, 
and  K2)  and  the  four  largest  diurnal  tidal  constituents  (Kl, 
01,  PI,  and  Ql).  We  use  the  parameterised  topographic 
internal  wave  drag  scheme  of  Garner  (2005),  modified  to 
limit  the  maximum  decay  rates  to  (9  hours)  ^  The  wave 
drag  parameterisation  represents  the  drag  resulting  from 
the  generation  and  subsequent  breaking  of  unresolved 
small-vertical  scale  internal  waves  by  tidal  flow  over  rough 
topography.  Because  the  wave  drags  acting  on  tidal  versus 
non-tidal  motions  have  different  strengths  (Bell,  1975),  tidal 
and  non-tidal  bottom  flows  are  separated  using  a  25-hour 
running  boxcar  filter.  Presently,  the  scalar  approximation 
(Ray,  1998)  is  used  for  the  self-attraction  and  loading  term 
(Hendershott,  1972). 

Barotropic  tidal  sea  surface  elevations  in  HYCOM  have  been 
compared  to  measurements  made  from  tide  gauges  (Shum 
et  al.,  1997)  and  state-of-the-art  data-assimilative  barotropic 
tide  models  (Egbert  et  al.,  1994).  In  Shriver  et  al.  (2012)  we 
also  compared  the  modelled  internal  tide  perturbation  to 
sea  surface  height  (SSH)  -  computed  from  a  spatial  high- 
pass  of  the  total  amplitude  of  tidal  SSH  -  to  internal  tide 
perturbations  computed  from  along-track  satellite  altimeter 


Figure  1.  The  internal  tide  amplitude  of  the  principal  lunar 
semidiurnal  constituent  M2  from  the  (a)  altimetric-based 
and  (b)  HYCOM  tidal  analyses.  The  five  subregions  denoted 
by  black  boxes  in  (b)  are  used  to  compute  area-averaged 
amplitudes  in  Shriver  et  al.  (2012).  From  Shriver  et  al.  (2012). 


data  (Ray  and  Byrne,  2010).  Figure  1  shows  global  maps 
of  the  M2  internal  tide  amplitude  in  HYCOM  versus  along- 
track  altimeter  data.  The  hotspots  around  locations  such 
as  Hawai’i,  the  Aleutian  Islands,  the  Tuamotu  Archipelago 
in  the  tropical  central  South  Pacific,  and  Madagascar  match 
up  reasonably  well  in  the  two  maps.  Furthermore,  the  root- 
mean-square  (rms)  perturbation  magnitudes  over  the  five 
hotspot  regions  delineated  by  boxes  in  the  bottom  panel  of 
Figure  1  agree  with  the  altimeter  rms  values  to  within  about 
20%.  Comparisons  of  tidal  currents,  and  their  vertical 
structure,  in  HYCOM  versus  historical  moored  observations 
were  made  in  Timko  et  al.  (2012;  2013). 

We  have  used  the  HYCOM  tidal  simulations  to  distinguish 
between  tidal  and  non-tidal  contributions  to  the  wavenumber 
spectrum  of  SSH  (Richman  et  al.,  2012).  The  slope  of  the 
wavenumber  spectrum  is  of  great  theoretical  interest,  as 
it  contains  clues  about  the  dominant  dynamics  of  low- 
frequency  flows  (e.g.,  LeT raon  et  al.,  2008;  Xu  and  Fu,  2011). 
Because  of  the  10-day  repeat  track  time,  current  state-of- 
the-art  TOPEX/JASON  satellite  altimeters  alias  tides  into 
longer  periods  (Parke  et  al.,  1987).  In  contrast,  because  our 
model  output  is  written  at  hourly  intervals,  we  can  easily 
separate  low-  from  high-frequency  signals  in  HYCCM.  In 
regions  such  as  the  Kuroshio,  where  low-frequency  motions 
are  more  energetic  than  tidal  motions,  the  wavenumber 
spectrum  of  SSH  is  dominated  by  low-frequency  motions. 
However,  in  regions  where  internal  tides  are  energetic,  for 
instance  near  Hawai’i,  high-frequency  motions  dominate 
the  high-wavenumber  end  of  the  wavenumber  spectrum. 

The  results  of  Richman  et  al.  (2012)  imply  that  internal  tides 
will  have  to  be  accurately  removed  from  data  taken  by  the 
planned  high-resolution  wide-swath  satellite  altimeter  (Fu 
et  al.,  2012)  before  low-frequency  oceanic  motions  can  be 
investigated.  The  accuracy  with  which  internal  tides  can 
be  removed  from  altimeter  signals  depends  on  the  degree 
of  non-stationarity  of  the  internal  tides.  Internal  tide  non- 
stationarity  in  HYCOM  is  investigated  in  Shriver  et  al.  (2014). 


3.  Insertion  of  topographic  wave  drag  into 
eddying  simulations 

In  Trossman  et  al.  (2013)  we  again  use  the  topographic 
wave  drag  parameterisation  of  Garner  (2005),  this  time 
applied  to  low-frequency  motions.  The  Trossman  et  al. 

(2013)  simulations  did  not  include  tides.  For  simplicity,  we 
refer  to  the  Garner  (2005)  scheme  here  as  a  “wave  drag” 
scheme  even  though  it  includes  effects  of  topographic 
blocking,  which  yields  low-level  turbulence,  as  well  as  lee 
wave  production  and  breaking.  Maps  of  energy  dissipation 
by  quadratic  bottom  drag  and  by  the  parameterised 
topographic  lee  wave  drag,  both  averaged  over  one  year, 
are  shown  in  Figure  2.  As  in  the  “offline”  estimates  of 
Nikurashin  and  Ferrari  (2011)  and  Scott  et  al.  (2011),  as 
well  as  the  offline  estimates  done  in  Trossman  et  al.  (2013), 
the  inline  estimates  shown  in  the  bottom  panel  of  Figure  2 
indicate  substantial  lee  wave  drag  dissipation  in  the  Southern 
Ocean.  The  globally  integrated  internal  lee  wave  drag 
dissipation  in  the  inline  Trossman  et  al.  (2013)  estimates 
is  0.4  TW,  comparable  to  that  seen  in  the  offline  estimates 
of  Nikurashin  and  Ferrari  (2011),  Scott  et  al.  (2011),  and 
Trossman  et  al.  (2013).  However,  as  anticipated,  there  is 
indeed  a  strong  feedback  between  the  internal  lee  wave 
drag  and  the  model  flows.  Adding  wave  drag  to  HYCOM 
substantially  changes  the  modelled  near-bottom  eddy  kinetic 
energy  and  stratification  fields  (T rossman  et  al.,  2013). 
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a)  Bottom  drag  [log^p(Wm  ^)] 
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Figure  2.  LoglO  of  the  bottom  and  wave  drag  terms  in  the 
energy  budget  of  1/12  eddying  HYCOM  simulation  with  ,  where 
each  drag  term  has  units  [W  m-2]:  (a)  quadratic  bottom 
boundary  layer  drag  and  (c)  parameterised  internal  lee  wave 
drag.  Adapted  from  Trossman  et  al.  (2013). 


4.  Ongoing  and  planned  work 

In  ongoing  and  planned  work  with  the  HYCOM  tide 
simulations,  we  are  comparing  the  SSH  frequency  spectra 
in  HYCOM  versus  tide  gauges,  as  a  measure  of  the  accuracy 
of  the  partition  of  modelled  low-  versus  high-frequency 
energy  (Savage  et  al.,  paper  in  revision).  We  are  estimating 
tidal  aliasing  errors  in  altimetrically-derived  low-frequency 
wavenumber  spectra  and  spectral  fluxes  (Savage  et  al., 
paper  in  preparation).  We  are  also  comparing  the  kinetic 
energy  frequency  spectra  in  HYCOM  versus  moored  current 
meter  records  (Muller  et  al.,  paper  in  preparation).  We 
are  comparing  both  low-  and  high-frequency  temperature 
variances  (which  are  related  to  available  potential  energies) 
to  temperature  variances  in  moored  historical  records 
(Luecke  et  al.  and  Bassette  et  al.,  respective  papers  in 
preparation).  To  further  elucidate  the  energetics  of  internal 
tides,  we  are  preparing  global  maps  of  the  baroclinic  tidal 
energy  fluxes  and  barotropic-to-baroclinic  tidal  energy 
conversions  (Ansong  et  al.  and  Buijsman  et  al.,  papers 
in  preparation).  The  latter  work  is  part  of  the  National 
Science  Foundation-funded  Climate  Process  Team  project 
“Collaborative  Research:  Representing  internal-wave  driven 
mixing  in  global  ocean  models”,  which  focuses  on  improving 
estimates  of  internal-wave  mediated  mixing  in  the  ocean, 
and  is  led  by  Professor  Jennifer  MacKinnon  of  the  Scripps 
Institution  of  Oceanography.  Finally,  our  long-term  goals 
for  the  HYCOM  tides  work  include  improvement  of  tidal 
accuracy  through  the  inclusion  of  data  assimilation  and 
better  estimates  of  the  self-attraction  and  loading  term,  and 
usage  of  the  global  HYCOM  tidal  solution  to  force  regional 
and  coastal  models  at  their  open-water  boundaries. 


In  ongoing  and  planned  work  on  parameterised  topographic 
lee  wave  drag,  we  will  investigate  whether  the  inclusion  of 
wave  drag  into  eddying  models  improves  the  comparison 
of  the  models  with  observations.  We  will  also  compare  the 
dissipation  predictions  of  the  Garner  (2005)  and  Bell  (1975) 
schemes  with  dissipation  inferred  from  microstructure 
observations.  Finally,  we  will  investigate  the  impact  of  a 
more  complete  parameterisation  of  the  vertical  deposition 
of  internal  lee  wave  drag  on  the  energetics  and  abyssal 
circulation  in  eddying  models.  These  projects  are  all  led  by 
co-author  Trossman  (papers  in  preparation). 
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